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Abstract 

Combining time-dependent density functional calculations for electrons with molecular dynamics 
simulations for ions, we investigate the dynamics of excited carriers in a (3, 3) carbon nanotube 
at different temperatures. Following an hv = 6.8 eV photoexcitation, the carrier decay is initially 
dominated by efficient electron-electron scattering. At room temperature, the excitation gap is 
reduced to nearly half its initial value after ~230 fs, where coupling to phonons starts dominating 
the decay. We show that the onset point and damping rate in the phonon regime change with 
initial ion velocities, a manifestation of temperature dependent electron-phonon coupling. 

PACS numbers: 81.07.De, 82.53.Mj, 73.22.-f, 73.63.Fg 
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Understanding the microscopic decay mechanism of electronic excitations is a challenging 
problem in nanotechnology. The fundamental hurdle is the requirement of unprecedented 
computational resources to treat the time evolution of electronic and ionic degrees of freedom 
in real time, from first principles, for physically relevant time periods. Carbon nanotubes 
are the ideal system to study these effects due to their well-defined structure, intriguing 
electronic properties and a high potential for application in future electronic devices P|. 
Additional interest has recently been triggered by observing optical emission induced by an 
electric current [^, radio- wave emission from aligned nanotubes working as an antenna Q], 
a„d potential apphcat.on of nanotubes in a .node-locked la.er g Whether studying the 
nature of excitonic effects affecting photoabsorption |6|], or attempting to identify the maxi- 
mum switching frequency of future field effect transistors or lasers, the key to progress is to 
understand the lifetime of excited carriers in carbon nanotubes. 

Time-resolved femtosecond pump-probe spectroscopy, applied to nanotube samples, has 
identified two regimes in the decay of carrier excitations [if . An ultra-fast decay channel has 
been associated with electron-electron scattering, which causes an internal thermalization 
of the electronic system within 200 fs. A slower decay channel has been associated with the 
excitation of ionic degrees of freedom at later times. Extracting information about carrier 
decay in specific, isolated nanotubes from those pump-probe data is not straight-forward, 
since samples contain a mixture of metallic and semiconducting nanotubes with different 
diameters and chiralities, all of which interact with the substrate. 

Here we present the first ab initio real-time calculation of the decay mechanism of elec- 
tronic excitations in narrow carbon nanotubes, which exhibit strong electron-phonon cou- 
pl,ng B a. ou. app..oaeh. we t.eat the t,.e dependence of th^elee.on.c degrees of 
freedom using time-dependent density 



parameterized TDDFT calculations 



'unctional theory (TDDFT) Unlike in previous 
we determine concurrently the ionic motion in the 
evolving charge density distribution by direct integration using molecular dynamics (MD) 
simulations within the force field given by the density functional theory (DFT). Our results 
provide not only an unbiased interpretation of available experimental data, but also identify 
the time scales associated with the dominant carrier decay channels and their temperature 
dependence. Our conclusion that the excitation gap is reduced by electron-electron scat- 
tering first, followed by electron-phonon interaction, is general and not restricted to narrow 
carbon nanotubes. 
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Alternatively, selected aspects of excited carrier dynamics could be described by summing 
up all diagrams corresponding to the electron-electron and electron-phonon many-body in- 
teractions. Under stead y-st ate conditions, where the Fermi golden rule and Boltzmann 
transport picture apply [ij], electron-phonon coupling can be determined by integrating 
the scattering matrix elements for transitions involving particular phonon modes. The cor- 
responding scattering time constant, including its temperature dependence, could then be 
determined using the inverse momentum-space integrals of these matrix elements . This 
steady-state approach, however, is not applicable in the sub-picosecond regime of interest 
here, since electron-electron scattering processes have been averaged out. Dynamics of elec- 
tronic decay channels only, ignoring coupling to phonons, could also be inferred from the 
imaginary part of the electron self-energy, determined using the GW approximation jl^ . 

Unlike these approaches, our calculations offer an unbiased insight into the short-time 
dynamics of electronic excitations, since electron and phonon scattering are treated on the 
same footing in real time within TDDFT-MD, using the FPSEID (First Principles Simulation 
tool for Electron Ion Dynamics) code ^|. Since the systems of interest are well above zero 
temperature, the ionic motion is described classically, with the forces acting on the ions 
given by the density functional theory in the local density approximation (LDA) [l6| . The 
electron- ion interaction is described using norm-conserving soft pseudopotentials [l7|, and 
the valence wave functions are expanded in a plane wave basis with a kinetic energy cutoff 
of 40 Ry. 

To allow for a realistic representation of phonons in isolated (3, 3) nanotubes, we arrange 
96-atom supercells, which extend over eight primitive unit cells of the tube, in a hexagonal 
geometry with an inter-wall separation of 5 A. We find that F point sampling of the small 
Brillouin zone is sufficient to describe the electron-electron and electron-phonon coupling 
well. To describe the dynamics at a given temperature, we randomize the initial ion velocities 
according to the Maxwell-Boltzmann distribution for that thermodynamical ensemble. The 
excited-state dynamics then proceeds in three steps. First, we simulate the initial optical 
excitation by promoting a valence electron to an unoccupied Kohn-Sham (KS) state. Next, 
by fixing this electronic configuration, we use constrained DFT to determine the electronic 
structure in that excited state. Finally, we continue a full, unconstrained TDDFT-MD 
simulation starting with these initial ion velocities and charge distribution. 

Figure Q] depicts the electronic structure and photoabsorption spectrum in the narrow 
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(3,3) carbon nanotube [la]. The electronic band structure of the primitive unit cell, shown 
in Fig.ma), agrees with published data |20| also in displaying the curvature effect, which 
modifies the Fermi momentum from the zone folding value at 2/3 (F — X). The photoab- 
sorption spectrum in Fig. Wljo) displays the response to light with the E-field polarized along 
the tube axis. Since depolarization effects normal to the tube are not important due to the 
small cross-section (2^, we used Fermi's golden rule directly to calculate the optical dipole 
matrix elements between occupied and empty KS eigenstates |2l| . 

We found two relevant photoabsorption processes, labelled "A" and "B" in Fig. [T] The 

ri 

absorption spectrum is dominated by a peak labeled "A" |19], which stems from many 
independent electron-hole transitions in a narrow energy range. A complex competition 
between the electronic decay channels, each with a different time constant, and phonon decay 
channels makes a clear distinction between electronic and phonon contributions problematic 
and will not be attempted here. In the following, we rather focus on the transition labelled 
"B" in Fig. ^ which, in spite of its smaller oscillator strength, is better defined due to the 
absence of nearby transitions with a large oscillator strength and allows a clear distinction 
between electronic and ionic decay channels. Our findings in terms of dominant decay 
channels, we believe, are of general nature, applicable to other photoabsorption processes in 
related systems. 

The time evolution of the electronic spectrum following a "B" excitation at t = is 
shown in Fig. |2l The time evolution of the KS eigenvalues indicates a rapid reduction of 
the electron-hole excitation gap j^l to less than half its initial value of 6.8 eV within 500 fs 
from the photoexcitation. The individual electron (hole) eigenvalues show an initial tendency 
to reach the bottom (top) of the conduction (valence) band very rapidly, manifesting the 
efficient non-radiative decay of the optical excitation towards the lowest lying electronic 
excitation. Under room temperature conditions, depicted in Fig. |21 and achieved by initially 
thermalizing the ionic motion at T = 300 K, we find a noticeable reduction of the decay 
rate ~230 fs after the photoexcitation. Subsequent decay may involve less efficient radiative 
transitions or energy transfer to phonons. 

We believe that the predicted reduction of the electron-hole energy gap could in principle 
be observed in time-resolved experiments. Light emission as a result of electron-hole recom- 
bination could generally be expected, since the computed optical matrix elements between 
the excited electron and hole states fiuctuate for a long time around a finite value. Even 
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though non-radiative decay of the excitation towards the lowest excitonic state is typically 
faster than radiative decay, emission could possibly be resonantly enhanced by the probe 
laser. 

To gain additional insight into the electronic decay, including possible radiative transi- 
tions, we investigated the time evolution of the matrix containing electron-hole dipole matrix 
elements. Immediately after the photoexcitation, the only non-vanishing matrix elements 
involve a field along the tube axis, same as in the initial photo-absorption. Later on, as 
phonons weakened the dipole selection rules, we observed matrix elements involving a field 
normal to the tube axis to increase, and eventually to reach similar values as those involving 
a field along the tube axis. As the electron-hole dipole matrix elements become isotropic in 
time, also the radiative decay should not show any preferential polarization. 

In spite of substantial changes in its electronic structure, the (3, 3) nanotube retains its 
cylindrical shape throughout the simulation. Close inspection of the vibrational spectra in- 
dicates enhancement of particular phonon modes following the photoexcitation, in particular 
at u = 556 cm~^, close to the ujrbm = 536 cm~^ value reported for the the radial breathing 
mode (RBM) of the (3,3) carbon nanotube ^j. Our conclusion that the RBM may be ex- 
cited by photoabsorption is also supported by recent femtosecond pump-probe experiments, 
which report exciting coherent phonons, including the RBM, in semiconducting single wall 
carbon nanotubes using resonant sub-10 fs pulses in the visible range 2^. 

The time-dependence of the photoexcitation gap, discussed in conjunction with Fig. |2l 
tells little about the nature and temperature dependence of the processes controlling the 
hot carrier decay. To obtain this information, we investigated the photo-decay in nanotubes 
with initial ionic temperatures of T = 77 K, 150 K, and 300 K. Witnessing to the numerical 
long-time stability of our TDDFT-MD simulations, we find the total energy to be conserved 
with a maximum deviation of about 10~^ eV/atom over the 0.5 ps long simulations. FigureEl 
depicts the key result, namely the time dependence of the potential energy of the nanotube 
under different initial conditions. 

Due to the finite size of our unit cell, the potential energy fluctuates in time, with larger 
fluctuations occurring at higher temperatures. Since ensemble averaging over may runs 
would necessitate unreasonably large computer resources, we invoked ergodicity and time 
averaged our potential energy profiles over 50 fs long time periods. Our results, displayed by 
the dotted lines in Fig. IHl indicate that independent of the initial temperature, the potential 
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energy decreases in time. After reaching a cross-over point, an initially gradual decrease 
becomes steeper. We find the decay rate to be nearly linear in both regimes, as shown by 
the dashed lines in Fig. |21 Whereas the decay rate in the early regime is rather independent 
of the initial temperature, its length decreases from ~360 fs at 77 K to ~230 fs at 300 K. 
Also, in the regime beyond the cross-over point, the decay rate decreases with increasing 
temperature, suggesting a slower energy flow from the electronic system to the ionic lattice. 
This is a clear indication of a temperature-dependent electron-phonon coupling. 

In contrast to the time dependence of the potential energy discussed above, the electron- 
hole energy gap, shown in Fig. |21 first shows a rapid, then a more gradual decrease. The 
cross-over points between the two decay regimes at 300 K, based on the electronic spectra 
in Fig. 121 and the potential energy in Fig. [HI nearly coincide at t~230 fs, suggesting the 
same underlying mechanism. Our finding that the initial decay is fast in Fig. |21 yet slow 
and temperature independent in Fig. El suggests that early decay regime is dominated by 
relaxations within the electronic system only. The opposite behavior is observed in the 
second regime, suggesting involvement of phonons. In view of the stronger electron-phonon 
coupling in the narrow (3, 3) nanotube, our estimated cross-over point between the regimes 
should still remain valid as a lower limit for wider nanotubes. In view of uncertainties 
associated with samples containing mixtures of nanotubes, the value t~200 fs, observed ^ 
in single- wall nanotubes with diameters close to 1.2 nm, lies very close to our prediction. 

Even though the time evolution of the charge densit y, u nderlving the density functional 
force field on the ions, is based on Ehrenfest dynamics [2J], our calculated forces do reflect 
the electron and ion dynamics at particular temperatures without relying on the adiabatic 
potential energy surface for a particular excited state. Deviations from Ehrenfest dynamics 
only occur, when the system cannot be associated with the adiabatic potential energy sur- 
face of an individual excited state. This would give rise to non- vanishing off-diagonal matrix 
elements of the time-dependent Kohn-Sham matrix, evaluated in the basis of the correspond- 
ing time-dependent eigenstates. We have indeed observed such non-vanishing off-diagonal 
matrix elements, albeit very small in magnitude, and only during short time periods in our 
simulation Q]. As another point worth mentioning, the electron-hole attraction in exci- 
tonic states, which is not reproduced by the DFT, is known not to change the total charge 
density much and thus not to affect the electron and ion dynamics. Thus, in stark contrast 
to the strong influence of ions on the excited carriers, which we consider explicitly, we find 
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the influence of electronic excitations on the ion dynamics to be neghgible. All the above 
arguments, along with the high degree of total energy conservation in our system discussed 
earlier, justify our use of Ehrenfest dynamics. 

In conclusion, our large-scale ab initio simulations of photophysical processes in carbon 
nanotubes suggest that the decay of photoexcited carriers occurs on a sub-picosecond time 
scale. An initial rapid nonradiative electronic decay is followed by a temperature-dependent, 
slower decay involving phonons. The duration of the electronic decay regime and the de- 
cay rate in the phonon-dominated regime decrease with increasing temperature, suggesting 
temperature-dependent electron-phonon coupling. 
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FIG. 1: (a) LDA band structure and (b) optical absorption spectrum of an isolated (3, 3) nanotube 
for a dipolc field along the tube axis. Two dipole-allowed transitions are labelled "A" and "B". 
The Fermi level lies at eV. The absorption spectrum (b) has been convoluted with a 0.3 eV wide 
Lorentzian. 
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FIG. 2: (Color online) Dynamics of hot-carrier decay following the excitation labelled "B" in 
Fig. after the ions have been thermalized at T = 300 K. Valence bands are shown by solid lines 
and conduction bands by dotted lines. The excited electron and hole states are shown by heavy 
solid lines. Dashed lines are guides to the eye. 
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FIG. 3: (Color online) Time dependence of the ionic potential energy after the photoexcitation 
at t = 0. The different solid lines describe systems with an initial Maxwell-Boltzmann velocity 
distribution corresponding to T = 77 K, 150 K, and 300 K. The light dotted lines represent 
time-averaged values. The dashed lines, showing linear fits to the data, are guides to the eye to 
distinguish the two time regimes. 
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